Transverse electric plasmons in bilayer graphene 
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We predict the existence of transverse electric (TE) plasmons in bilayer graphene. We find that 
their plasmonic properties are much more pronounced in bilayer than in monolayer graphene, in a 
sense that they can get more localized at frequencies just below fkJ = 0.4 eV for adequate doping 
values. This is a consequence of the perfectly nested bands in bilayer graphene which are separated 
O ' by ~ 0.4 eV. 
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' Plasmons arc self-sustained collective electron excitations which are of great interest both for fundamental physics 
r T | \ and potential technological applications. Plasmon is a paradigmatic quantum many-body phenomenon studied in 
condensed matter physics [l|. Closely related excitations are surface plasmons which hold promise as a possible tool 
for controlling light at subwavelength scales 041), giving rise to the field of nanophotonics, and they also play an 
important role in metamaterials |6H9j; plasmons in essentially two-dimensional (2D) structures are similar in this 
respect to surface plasmons. These reasons are a great motivation for studying plasmonic excitations and their 
properties in novel materials. Two such materials are monolayer [l(| EH and bilayer (see e.g., [ill, E3) graphene. 
Graphene is a 2D sheet made of carbon atoms organized in a honeycomb lattice [lfj . Ill| , whereas bilayer graphene 
pi I consists of two such layers stacked on top of each other in a certain way [ill, E3 ■ While there are only a few studies 
O ■ of plasmons in bilayer graphene [13l4l6| . these collective excitations have attracted substantially more attention 
^ | in monolayer graphene [1714231 ]. Several years ago it was predicted that graphene, besides the ordinary longitudinal 
plasmons [transverse magnetic (TM) modes] fl7H20l . [23| . also supports unusual transverse plasmons [transverse electric 
' S ' (TE) mode s] | 19| . These excitation are possible only if the imaginary part of the conductivity of a thin sheet of material 
is negative [l9( • On the other hand, such a conductivity requires some complexity of the band structure of the material 
i-C ' involved. For example, TE plasmons cannot occur if the 2D material possesses a single parabolic electron band. From 
this perspective, bilayer graphene, with its rich band structure and optical conductivity (e.g., see [24[ and references 
therein), seems as a promising material for exploring the possibility of existence of TE plasmons. Here we predict the 
t-H | existence of TE plasmons in bilayer graphene. We find that their plasmonic properties are much more pronounced in 
^ i bilayer than in monolayer graphene, in a sense that the wavelength of TE plasmons in bilayer can be smaller than in 
CD ' monolayer graphene at the same frequency. 

Throughout this work we consider bilayer graphene as an infinitely thin sheet of material with conductivity er(q, lu). 
We assume that air with e r — 1 is above and below bilayer graphene. Given the conductivity, by employing classical 
electrodynamics, one finds that self-sustained oscillations of the charge occur when (see [i~9j and references therein) 
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for TM modes, and 

1 _ /i a;zcr(q,u) = Q ,^ 
H : 2vV - ^ 2 /c 2 
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for TE modes. The TM plasmons can considerably depart from the light line, that is, their wavelength can be 
considerably smaller than that of light at the same frequency. For this reason, when calculating TM plasmons it is 
desirable to know the conductivity as a function of both frequency w and wavevector q. However, it turns out that 
the TE plasmons (both in monolayer [T9[ and bilayer graphene, as will be shown below) are quite close to the light 
line q = uj/c, and therefore it is a good approximation to use <j(u>) = cr(q = 0,u;). Moreover, these plasmons are 
expected to show strong polariton character, i.e., creation of hybrid plasmon-photon excitations. At this point it is 
worthy to note that if the relative permittivity of dielectrics above and below graphene are sufficiently different, so 
that light lines differ substantially, then TE plasmon will not exist (perhaps they could exist as leaky modes). 

The conductivity a(uj) = $ia(ui) + i^sa(cu) is complex, and plasmon dispersion is characterized by the imaginary 
part Skt(u;), whereas 5Rcr(w) determines plasmon losses, or more generally absorption of the sheet. From Eq. ([2]) it 
follows that the TE plasmons exist only if ^sa(uj) < |19j . 

In order to calculate the imaginary part of the conductivity, we employ Kramers-Kronig relations and the calculation 



of absorption by Nicol and Carbotte [24|, where ?R.a(oj) [see Eqs. (19)-(21) in Ref. [24[] was calculated by using the 
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FIG. 1: The band-structure of bilayer graphene. The two upper bands (as well as the two lower bands) are perfectly nested 
and separated by 7 ~ 0.4 eV; qo = j/hvF- Horizontal line depicts one possible value of the Fermi level, and arrows denote 
some of the possible interband electronic transitions. See text for details. 



Kubo formula. The optical conductivit 
is organized in four bands given by 



has rich structure due to the fact that the single-particle spectrum of graphene 
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where vp = 10 6 m/s, the parameter 7 0.4 eV is equal to the separation between the two conduction bands (which 
is equal to the separation between the valence bands). The band structure © is calculated from the tight binding- 
approach, where vf is connected to the nearest-neighbour hopping terms for electrons to move in each of the two 
graphene planes, and the distance between Carbon atoms in one monolayer (see Ref. 24 Q, whereas 7 is the hopping 



parameter corresponding to electrons hoping from one layer to the other and vice versa [241 ] . The two graphene layers 
are stacked one above the other according to the so-called Bernal-type stacking (e.g., see Ref. (HI). We emphasize 
that the perfect nesting of bands gives rise to the stronger plasmon like features of TE plasmons in bilayer than in 
monolayer graphene. The four bands are illustrated in Fig. [T] along with some of the electronic transitions which 
result in absorption. Absorption depends on 7 and the Fermi level /z; the latter can be changed by applying external 
bias voltage. 

The imaginary part of the conductivity can be calculated from 3ia(uj) by using the Kramers-Kronig relations 



which yields 
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an = e 2 /2h, 0(x) = 1 if x > and zero otherwise, and f2 = hui. Here we assume zero temperature T ~ 0, which 
is a good approximation for sufficiently doped bilayer graphene where \i ^> fcsT. Formulae ([5]) and (|6]) are used to 
describe the properties of TE plasmons. 

In Figure [5] we show the real and imaginary part of the conductivity for two different values of the Fermi level: 
fj, = O.47 and /1 = O.97 (wc focus on the electron doped system fi > 0). Because plasmons are strongly damped 
by interband transitions, it is instructive at this point to discuss the kinematical requirements for the excitation of 
electron- hole pairs. If the doping is such that /i < 7/2, a quantum of energy (plasmon or photon) with in-plane 
momentum q = can excite an electron- hole pair only if Hlj > 2\i (excitations from the upper valence to the lower 
conduction band shown as red dot-dashed line in Fig. [T]). If \i > 7/2, the (q = 0, cj)-quantum can excite an electron- 
hole pair only for huj > 7 (excitations from the lower to the upper conduction band shown as green solid lines in Fig. 
Q] occur at Huj = 7). If the plasmon/photon has in-plane momentum q larger than zero, then interband transitions 
are possible for smaller frequencies (see blue dashed lines in Fig. [IJ. There is a region in the (g,w)-plane where 
electron- hole excitations are forbidden due to the Pauli principle (e.g., see figures in Refs. fl4T - [il| ). Because plasmons 
are strongly damped by these interband transitions (this is Landau damping), in our search for the TE plasmons, 
we focus on their dispersion curve in the regime where electron-hole pair formation is inadmissible (via first-order 
transition). 




FIG. 2: The real (red dotted lines) and imaginary (blue solid lines) part of the conductivity of bilayer graphene for two values 
of doping: /i = O.47 (a), and fi = O.97 (b). The conductivity is in units of ao = e 2 /2h, and the frequency is in units of cuo = f/H. 
The 5-functions in 5Ra(oj) at ui — (intraband transitions) and u) = 7/fi, (transitions from the lower to the upper conduction 
band depicted as green solid arrows in Fig. [1]) are not shown (see [241]). 

In Figure [3] we show the plasmon dispersion curves for fi = OAj and /i = O.97; in the spirit of Ref. [T^], we show 
Aq = q — uj/c as a function of frequency u). Plasmons are very close to the light line and thus one can to a very good 
approximation write the dispersion curve as 

A^^(^. (7) 

To the left (right) of the vertical red dotted line in Fig. |3l plasmon damping via excitation of electron-hole pairs is 
(is not) forbidden. 
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FIG. 3: The plasmon dispersion curve Aq = q — ui/c vs. uj for ^ = O.47 (a), and \x = O.97 (b) is shown as blue solid line. 
To the right of the vertical red dotted lines plasmons can be damped via excitation of electron-hole pairs, whereas to the left 
of this line these excitations are forbidden due to the Pauli principle. Black dashed line in (b) (which closely follows the blue 
line) corresponds to Eq. JHJ. The wave vector is in units of qo = y/hvF, and the frequency is in units of uio = j/H. 



For [i = O.47, 3<t(w) is smaller than zero for ui in an interval of frequencies just below 2/z. From the leading term 
in 3cr(u;) we find that departure of the dispersion curve from the light line is logarithmically slow: A<7 <^< 7 /2 oc 
[log \Hui — 2//|] 2 . The same type of behavior occurs in monolayer graphene fl9| . 

However, for fj. = O.97, one can see the advantage of bilayer over monolayer graphene in the context of TE plasmons. 
The conductivity ^sa(ui) is smaller than zero in an interval of frequencies below 7. In this interval, the most dominant 
term to the conductivity is the last one from Eq. ([5]). that is, 

A<7 7/2<Al < 7 « 

This approximation is illustrated with black dashed line in Fig. [31 and it almost perfectly matches the dispersion 
curve. Note that the singularity in ^sa(ui) at tuo = 7 is of the form 1/(7 — Tia;), whereas the singularity at frio = 2/j, is 
logarithmic (as in monolayer graphene [ijj ). As a consequence, the departure of the dispersion curve from the light 
line in bilayer graphene is much faster for /1 > 7/2 than for /1 < 7/2, and it is faster than in monolayer graphene 
as well [note the two orders of magnitude difference between the abscissa scales in Figs. [21a) and (b)]. Thus, we 
conclude that more pronounced plasmonic features of TE plasmons (shrinking of wave length which is measured as 
departure of q from the light line) can be obtained in bilayer graphene. The term in 3yt(w) which is responsible for TE 
plasmons for \x > 7/2 corresponds (via Kramers-Kronig relations) to the absorption term b(/j,)S(hui — 7) (24|, which 
arises from the transitions from the first to the second valence band (shown as green solid arrows in Fig. [1]), which 
are perfectly nested and separated by 7. Thus, this unique feature of bilayer graphene gives rise to TE plasmons with 
more pronounced plasmon like features than in monolayer graphene. 

Before closing, let us discuss some properties and possible observation of TE plasmons. First, note that since the 
electric field oscillations arc both perpendicular to the propagation vector q, and lie in the bilayer graphene plane, the 
electric current j = er(u;)E is also perpendicular to q. Thus, j • q = 0, and the equation of continuity yields that the 
charge density is zero (i.e., one has self-sustained oscillations of the current). In order to excite plasmons of frequency 
ui with light of the same frequency, one has to somehow account for the conservation of the momentum which is larger 
for plasmons. Since the momentum mismatch is relatively small, the standard plasmon excitation schemes such as 
the prism or grating coupling methods (e.g., see and references therein) could be used for the excitation of these 
plasmons. 

In conclusion, we have predicted the existence of transverse electric (TE) plasmons in bilayer graphene. Since they 
exist very close to the light line, these plasmons are expected to show strong polariton character, i.e., mixing with 
photon modes. However, due to the perfectly nested valence bands of bilayer graphene, their dispersion departs much 
more from the light line than in monolayer graphene. 
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